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Recent experiments on high density excitons in CU2O provide evidence for degenerate quantum 
statistics and Bose-Einstein condensation of this nearly ideal gas. We model the time dependence 
of this bosonic system including exciton decay mechanisms, energy exchange with phonons, and 
interconversion between ortho (triplet-state) and para (singlet-state) excitons, using parameters for 



^\ , the excitonic decay, the coupling to acoustic and low-lying optical phonons, Auger recombination, 

OS ■ 

and ortho-para interconversion derived from experiment. The single adjustable parameter in our 
model is the optical-phonon cooling rate for Auger and laser-produced hot excitons. We show that 
the orthoexcitons move along the phase boundary without crossing it (i.e., exhibit a "quantum 
saturation"), as a consequence of the balance of entropy changes due to cooling of excitons by 
phonons and heating by the non-radiative Auger two-exciton recombination process. The Auger 
annihilation rate for para-para collisions is much smaller than that for ortho-para and ortho-ortho 
collisions, explaining why, under the given experimental conditions, the paraexcitons condense while 
the orthoexcitons fail to do so. 



PACS numbers: 25.70.Np,12.38.Qk 



I. INTRODUCTION 



While Bose-Einstein condensation underlies the remarkable properties of the strongly-interacting superfluid, liquid 
helium II, discovery of Bose-Einstein condensation in nearly ideal gases has proven difficult. Many recent exper- 
iments have focussed on this challenge In the semiconductor CU2O, excitons (electron-hole pairs bound by 
their Coulomb attraction) have exhibited Bose-Einstein statistics ||,[| and indeed Bose-Einstein condensation HQ. 
In these experiments, intense pulses of laser light excite the crystal, creating a gas of (spin-triplet) orthoexcitons 
I and (spin-singlet) paraexcitons, split by an exchange energy AE w 12 meV H||. The kinetic energy distribution of 
orthoexcitons as a function of time from the onset of the laser pulse is observed by spectroscopy of their photolumi- 
nescence (considerably more intense than that of the paraexcitons). In the classical (low density) regime, the energy 
distribution is observed to be Maxwell-Boltzmann, described by an effective exciton temperature. In the quantum 
(high density) regime, the spectra exhibit Bose-Einstein distributions, as shown in Figs. l(a)-(c). These spectra 
are well fit in terms of an ideal gas with an instantaneous chemical potential, /x, and temperature, T. From these 
parameters, one can directly calculate the density of the gas. 

The densities obtained from the recombination spectra of Snoke et al. are plotted in Fig. 1(d). This figure shows 
the experimental results with a long (10 ns) laser pulse (open circles) and with a short (100 ps) laser pulse (black 
dots). These experiments observe a "quantum saturation" of the (spin-triplet) orthoexcitons, i.e., a tendency for them 
to move closely parallel to the critical line without condensing The critical line, an adiabat (constant entropy per 
particle, s = S/N), has the form, T c = {2-rth 2 /mk B )(n/ g((3/2)) 2 / 3 , where the de generacy g= 3 for orthoexcitons and 
1 for paraexcitons. We note that lines in the phase diagram parallel (in a log-log plot) to the condensation line at 
higher temperature correspond to adiabats, along which a = —n/ksT is constant, since s is only a function of a for 
an ideal Bose gas. 

Experimental measurements of the paraexcitons indicate condensation of this component in unstressed || and 
uniaxially stressed crystals [Q. The stressed results are summarized in Figs. 2(a,b). At high densities the orthoexcitons 
decay very rapidly, due, as we shall see, to Auger recombination and ortho-to-para conversion. In contrast, the 
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paraexcitons decay very slowly even at high densities. Figure 2(b) shows density-temperature trajectories for both 
orthoexcitons (open circles) and paraexcitons (black dots) for long pulses. The stress splits the triply degenerate 
orthoexciton level into three components and only the lower of the three is significantly populated, leading to a closer 
proximity to the condensation line than in the unstressed case. No orthoexciton condensation is observed. The 
stress, however, does lead to clear evidence in the paraexciton spectrum that the paraexcitons condense: their energy 
distribution is much sharper than that of the orthoexcitons at a given time. Under these conditions the paraexciton 
density is determined from the relative intensities of ortho and paraexcitons, combined with the spectroscopically 
determined density of the orthoexcitons. This analysis yields the striking evidence for paraexciton condensation 
shown in Fig. 2(b). 

The non-stressed paraexciton spectrum also displays anomalies which are interpreted as a result of condensation 
but are not so cleanly analyzed: a highly blue-shifted component is observed which may indicate a superfluid with 
a large drift velocity. The exact effect of the stress is not well understood, but it is possible that it helps to localize 
better the two components of the gas, making the analysis more straightforward Jul. 

In this study we identify the salient factors that lead to the quantum saturation of orthoexcitons. We find, by 
quantitative investigation of the system, that the observed behavior can be understood in terms of well-established 
kinetic processes of cooling, reheating, ortho-para conversions, and recombination of the excitons. The essence of 
the argument is that the rate of change of the entropy of the excitons is a balance between the entropy loss due to 
phonon cooling and heating following Auger annihilation of exciton pairs into hot electrons and holes. The rate of 
phonon cooling, which is dominated by emission of acoustic phonons, varies as aT 3 / 2 for exciton temperatures large 
compared with the ambient lattice temperature. On the other hand heating of the excitons from the Auger process 
has, as we shall see, the form bn a proportional to the density of orthoexcitons, n a . Thus the time rate of change of 
the orthoexciton entropy per particle, s a , is given in first approximation by 

= -aT 3 / 3 + bn . (1) 

at 

Note that the right side vanishes along the adiabat given by T 3 / 2 /n D — b/a, which is a stable fixed point of this 
equation, e.g., if the temperature is too high, phonon cooling will increase, pulling the entropy down to this adiabat, 
etc. As we see in Sec. Vic, where we present this argument in greater detail, the stable adiabat is only slightly above 
the condensation line. Our qualitative and quantitatively detailed explanation of the quantum saturation indicates 
that we have identified the critical mechanisms of the exciton dynamics, enabling us to assert with confidence that 
the paraexcitons do indeed condense under the experimental conditions. 

We model the time evolution of the gas from a rate equation approach. In Sec. II we construct the rate equations 
that incorporate the relevant dynamical processes, which we examine individually in Sees. III-V. In Sec. Ill we consider 
exciton-phonon interactions and show how the onset of quantum statistics tends to suppress these cooling processes. 
The effects of ortho-para conversion mechanisms |n],|l2| are considered in Sec. IV. We find that the ortho-para 
conversion process is strongly modified for excitons in the condensed phase. We examine in Sec. V the Auger-decay 
process, in which two excitons annihilate to form an ionized electron-hole pair p3| , with the ionized electrons and holes 
subsequently recombining into excitons. This process tends to keep the gas away from the degenerate region both by 
decreasing the exciton density and by indirectly heating the excitons. At high densities the orthoexcitons decay very 
rapidly due to Auger recombination and ortho-to-para conversion; however, as we show, the striking longevity of the 
paraexcitons, which is crucial to their ability to cross the Bose-Einstein condensation boundary as seen in Fig. 2(b), is 
due to the weakness of the para-para Auger process. This weakness is a result of the detailed band structure. In Sec. 
VI an analysis of the solutions to the rate equations for a range of values of the input parameters allow us to see the 
relative importance of the various physical processes in producing the quantum saturation and condensation effects. 
In the same section we discuss the evolution of the orthoexciton entropy, and calculate the time for the orthoexcitons 
to approach the adiabat of stability. Finally in Sec. VII we summarize our conclusions. 
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II. RATE EQUATIONS 



We treat the exciton gas as nearly ideal, and kept in thermal equilibrium by the rapid exciton-exciton collisions, 
with a single temperature T, generally greater than the lattice temperature Tg. For example, for excitons scattering 
as hard spheres with a scattering radius on the order of the exciton Bohr radius, the characteristic scattering time is 
a few picoseconds for typical densities 10 18 -10 19 cm~ 3 . The slower interconversion processes between the ortho and 
paraexcitons do not allow them to come into chemical equilibrium quickly, so that, in general, the chemical potential 
of the orthoexcitons /x differs from that of the paraexcitons \i p . We further assume that the lattice maintains a 
constant temperature Tg_ throughout the relaxation process. At 30-50 ns after an intense 10 ns pulse the excitons 
assume a temperature ~ 5-6K, which we conclude is the lattice temperature in the excitation region at that time, 
somewhat higher than that of the surrounding cold lattice and the bath temperature (2K). 

In the experiments, photoexcitation by green laser light creates hot electrons and holes within about a micrometer 
of the crystal surface. These hot carriers diffuse and within a nanosecond relax to form the excitons under study. 
During the first 5-20 ns the spatial extent of the excitonic gas under consideration has been experimentally estimated 
to be about 30 /im. We make the simplifying assumption that the excitonic gas occupies a constant volume over 
the entire time interval, which is valid at least for the first 7 ns. Going beyond this approximation to include the 
full hydrodynamics of expansion is a task for a future paper (see jl4| for an initial approach). We also neglect 
inhomogeneities in the exciton gas and assume that the ortho and paraexcitons occupy the same volume. The 
parameters of the problem, all time-dependent, are therefore the exciton temperature T, the number of orthoexcitons 
iV , and the number of paraexcitons N p . 

In the normal (non-condensed) regime the chemical potentials fa (i = o,p) of each of the species are functions 
of Ni and T. In the condensed regime, where the chemical potential is zero, we take the number of condensed 
particles Ni, con d to be the equilibrium result, Ni tCon d = iV»(l — (T/T C j) 3 / 2 ), where T C i is the critical temperature for 
condensation of the component i. 

Excitons are formed by the laser through creation of ionized electrons and holes, with an excess kinetic energy 
of about 250 meV per pair, which then combine in random spin states to give orthoexcitons and paraexcitons in a 
three-to-one ratio. The rate equations that describe how the total numbers N a and N p of the two species change in 
time are then: 



= G„(t) -zr- DN , exc + UN p>exc - + + ) (2) 
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where Gi(t) is the laser production rate of excitons, and the Tu are the intrinsic radiative lifetimes (t; c sa 300 ns and 
n p « 150 lis). The quantity D is the ortho to para down-conversion rate, and U is the para to ortho up-conversion 
rate; as we show later only the non-condensed (non-zero momentum) excitons participate significantly in the ortho- 
para interconversion processes, so that these terms are effectively proportional to the number of non-condensed, or 
excited, particles, Ni^ exc . The terms —Ni/rAi describe the Auger annihilation process, while the final terms account 
for the reformation of excitons from the Auger- ionized electrons and holes. The factor of 1/2 arises because for each 
two excitons destroyed by the Auger process, only one reforms, with probability 1/4 to be a paraexciton and 3/4 to be 
an orthoexciton. The right side of these equations is a function of the numbers of excitons N a , N p and their common 
temperature T. 

The third kinetic equation is the rate of change of the total internal energy, E = E Q + E p , of the gas. 



dE\ _ /dE\ /dE\ /dE\ f dE\ f dE\ / dE 

"til ~ \^) ph onon + \^)o^p + \^)p^o + \^)L + \^)A + \M 



(4) 



The terms in this equation describe, respectively, the change in the energy due to collisions of the excitons with the cold 
lattice ("phonon"), phonon-assisted ortho-para interconversions, laser heating (L), heating from Auger annihilations 
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(A) and radiative recombination (1). In fact, because the radiative recombination lifetimes ti q and n p are long, the 
latter process is negligible in the three rate equations. Since E is a function of N Q , N p and T, the rate equations can 
be equivalently cast as equations for the rate of temperature change dT/dt, and the exciton numbers. Let us now 
examine each of the relaxation mechanisms in turn. 



III. PHONON COOLING OF EXCITONS 

Excitons in CU2O interact with acoustic and optical phonons through a deformation potential interaction |l^,|l6|. If 
the gas is above the condensation temperature T c , the rate of change of the energy of the exciton gas for the processes 
in which an exciton emits or absorbs a phonon, shown in Figs. 3(a,b), is: 

TiT ) =tI1 l M d 2 W(i + - (1 + »jg)(l + /.r)}^? - H - ^<f); (5) 

Ul / phonon r - 

here is the (ortho or para) exciton distribution, is the phonon distribution, is the exciton-phonon interaction 
matrix element, Sj: is the exciton energy and hwq the phonon energy. Below T c , the ground state is occupied by a 
macroscopic number of particles, so that we have an additional contribution, Figs. 3(c,d), from condensed particles: 

dE\ con< ^ 27r 

-XT ) = "? £ \M<j\ 2 {N cond (l + n ? )ff - (1 + N cond )n ? (1 + / r ) }fo«; r *(e ? - (6) 

^ / phonon 11 g 

We define the exciton gas structure factor, 

S(<f, w) = ^ ng(l + «fe +f )<5(£fc +f - £% - M, (7) 

k 

in terms of which Eq. (||) takes the form: 

§) = 2 i E |Md 2 ^(g> ? -)/,(l - e^>^), (8) 

OT / phonon 11 g 

where (3t = l/ksTg and f3 = l/fc^T. In terms of u± = f3[(h 2 /2m)(q/2 ± mu/hq) 2 — /j], the structure factor becomes 

V m 2 1 / 1 - e~ v + 



= — 71 — 73l7JT ln (l ^7-)' ( 9 ) 



where V is the volume of the exciton gas and m is the exciton mass, w (3.0 ± 0.2)m° (where m° is the bare electron 
mass) 0,|8| for the orthoexcitons; in this analysis we assume that both the ortho and the paraexcitons have the 
same massj^] and let eg = E g — Et, + AE + h 2 k 2 /2m, where E g is the energy gap energy (w 2.17 eV in CU2O), Eb 



x Our choice of the effective masses is based on the experimentally known Rydberg of 97 meV for the excited exciton states, 
which depends on the reduced mass of the electron and the hole and the dielectric constant of the material. A range of values 
for the effective electron and hole masses have been reported [po|]2l| l , and there are also uncertainties in the dielectric constant 
[[l9| . The present calculation is not critically sensitive to the choice of masses and dielectric constant, so we choose the measured 
values m e = 0.84mg, = 0.61m° and eo = 7.11 [[L9| for the static dielectric constant, which yield the excitonic Rydberg w 
99 meV, very close to the experimental value obtained from the absorption spectrum of the excited excitonic states. The total 
orthoexciton mass is not equal to the sum of the effective electron and hole masses rm, likely due to the central cell corrections 
, i.e., the corrections due to fact that the Bohr radius is comparable with the lattice constant. The central cell corrections 
are also likely the reason that the binding energy of the Is state (approximately 153 meV) is much larger than the Rydberg for 
the excited excitonic states. This measured binding energy corresponds to a Bohr radius a = e 2 /2eo-Et ft* 7A, where Eb is the 
exciton binding energy. 
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is the binding energy, and AE is the exchange interaction (in which an exciton virtually annihilates and reforms), 
which is non-zero for the orthoexcitons and zero for the paraexcitons. 

For acoustic phonons u>$ = v^q, where vi is the longitudinal sound velocity. We neglect the interaction of excitons 
with acoustic transverse phonons, since it is about fifty times weaker than that with the longitudinal phonons [ ]l7| , p3|] . 
The square of the matrix element for deformation potential scattering is (Tq | 



1 9 1 2pVv t (l + (ga/4) 2 ) 4 ' 



where D a is the acoustic deformation potential, « 1.8 eV experimentally p8[ , a is the exciton Bohr radius, q is the 
phonon wave vector, and p « 6.0 gm/cm 3 is the mass density of the material. The quantity (1 + (ga/4) 2 ) 4 in the 
denominator of Eq. (|l^) comes from the square of the Fourier transform of the relative Is electron-hole wave function, 
evaluated at m e = m^, where are the effective electron and hole masses. This correction term becomes significant 
when the wavelength of the phonon is comparable to the Bohr radius, and is actually negligible for thermalized 
excitons with T < 7 OK, cf. |f|. 

Combining Eqs. (||) and (fio|), we find the rate of acoustic-phonon cooling per exciton, 



1 (dE\ _ (mvj) 2 1 m 1 f a\ 3 {k B TfD 2 a 



N \ dt J LA fi 8ir 3 mi on na 3 \d£ J (hvi/ai) 6 (mv 2 ) 

1 1 



dx x 



1 e xT/T e _ Y 



(1 + {k B Tax/Ahv e ) 2 ) 4 V 1 - e ~ y - 



where n — N/V is the density of excitons, ai = 4.26A is the lattice constant, rrii on = paf is a mean ion mass and 
a w 7 A for the ground state. The factors v± in Eq. (g) become y± — (mv 2 /8kBT)[(xkBT/mv 2 ) ± 2] 2 — Pp. The 
prefactor (mv 2 ) 2 /h is 176 meV/ns. Equation ( pT| ) depends upon the quantum degeneracy factor, a = —p/ksT. 
Along lines of small constant a, the acoustic-phonon cooling rate per particle, with D a = 1.8 eV from Snoke et al., is 



1_ fdE_ 
N \~dt 



w -0.49T 3/2 (1 - Te/T) meV/ns, (12) 

LA 



with T measured in Kelvin; in the classical limit, the coefficient becomes 0.62. 

Excitons also cool by emission of optical phonons, of which there are 15 modes at the zone center of CU2O p2[] . The 
optical phonons are: r^" 5 (three-fold degenerate, of energy 11.4 meV), r^ 2 (two-fold degenerate, of energy 13.8 meV), 
r^ 5 (three-fold degenerate, of energy 18.7 meV), (non-degenerate, of energy 43 meV), r^~ 5 (three-fold degenerate, 
of energy 64 meV), and T^ 5 (three-fold degenerate, of energy 79 meV). They group into two parts; the first eight 
with relatively lower energies, while the remaining seven have relatively higher energies. In calculating the cooling 
rate of thermalized excitons we neglect the contribution of the second group; even the lowest ones with Tim = 43 
meV correspond to a Maxwellian temperature of 2huj/3kB ~ 300K, which suppresses their contribution to exciton 
cooling. In contrast, these high energy phonons definitely need to be taken into account in calculating cooling of non- 
thcrmalizcd high kinetic energy excitons, Sec. V. All the lower energy phonons have negative parity, as their group 
theory notation indicates. Since the parity of the excitons is even, the matrix element for the exciton odd-parity 
optical phonon interaction vanishes at the zone center. The average squared matrix element for the deformation 
potential interaction with an optical phonon of frequency u>i can be written as |]l6| 

Wf\ 2 = (13) 

where Di is the optical phonon deformation potential, averaged over the whole Brillouin zone; for the T^ 5 phonon of 
zone-center energy 18.7 meV, Di w 0.17 eV/A in magnitude [fl8|| . As noted in Ref. |l8) , the ratios of the deformation 
potentials of the low-energy phonons can be set by comparison with the data of Yu and Shen (lj] for the relative 
efficiency of the three "thrcc-phonon Raman lines." The heights of these lines, which are simply due to hot lumines- 
cence following emission of an optical phonon by an exciton, are proportional to the square of the nonpolar optical 
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deformation potentials involved, divided by the phonon energy. In the data of Ref. fl7| ], the r]" 5 phonon emission line, 
which corresponds to the sum of the contributions of two TO (transverse optical) and the LO (longitudinal optical) 
mode, has a maximum intensity about three times higher than the two lower-energy phonon lines, which are roughly 
equal in maximum efficiency. Using this analysis Snoke et al. have fixed the deformation potentials of the r^" 5 and 
the rj" 2 phonons at 0.72 ± 0.08 times the T^ 5 deformation potential. 

The contribution from the optical phonons to the cooling rate per exciton is, 

— (—\ \ " ( mv i) 2 1 m 1 (k B T) 2 (Dia) 2 ( 1 1 

~~N \ dt ) opt ~ ^ h tt3 m ion no? (hv e /ai) 4 \ e^ h ^ ~ 1 e ^ h ^ _ 

(M) 

where now the factors v± become w± = (x ± hu>i/AkBTx) 2 — f3fi. The sum is over the group of lower-energy optical 
phonons. Along lines of constant a, and small lattice temperature, the cooling rate per particle due to each of the 
optical phonons varies as j 11 /^-/ 3 ^/ 2 _ 

From Eq. (||), the net additional energy loss per condensed exciton, from excited excitons emitting an acoustic 
phonon and dropping into the condensate is 

1 fdE\ cond 8( m vff m (mvf)Dl 1 

{nk -Jk ), (15) 



N C ond \dt ) LA Trh m ion (hvi/ai) 3 (i + (k a/4) 2 ) 

where fco — 2mvi/h. For emission of optical phonons, 
1 

Ncond V dt J opt ^T 1 2 1/27r h rn lon \ a J (■mvj) 3 / 2 (hv e /ac 



dEV° nd ^ {mvjf m (D t a) 2 (h^/ 2 _ 



where k, t — (2muji/h) 1 / 2 . 

Combining these various rates, we find the rate of change of the internal energy of the excitons per particle due to 
the exciton-phonon interaction shown in Figs. 4(a,b) for acoustic and optical phonons, respectively. In these figures 
the temperatures of the exciton gas and the lattice are kept fixed at 30K and 6K, respectively, as the chemical 
potential and therefore the density is varied. As we see, for typical temperatures and densities, cooling of excitons 
by acoustic phonons dominates cooling by optical phonons by more than an order of magnitude. Optical phonons do 
not contribute significantly to the cooling since their energy is greater than that of the majority of the excitons. In 
contrast, optical phonons play the dominant role in Auger heating, as we discuss in Sec. V, since the process involves 
either free electrons and holes, or excitons of very large kinetic energies. 

The results in Figs. 4(a,b) show a clear effect of degeneracy on the exciton cooling. The rate is constant when 
the exciton gas is in the classical regime, but it decreases with increasing degeneracy, approaching a lower value as 
a — > 0. This decrease occurs because energy and momentum conservation limits the number of excitons that are able 
to interact with phonons as the degeneracy increases. As the gas approaches the condensation phase boundary, the 
cooling mechanisms become less effective. Furthermore the dependence of the acoustic phonon matrix element on 
momentum transfer (Eq. ([l0])) enhances this effect. As we discuss in the last section, the smaller cooling rate at high 
degeneracy plays a significant role in the evolution of the system. A related classical phenomenon has been observed 
by Trauernicht et al. p3( ] in experiments on CU2O, where the diffusion constant of paraexcitons shows a sharp rise at 
low temperatures due to their average velocity falling below the sound velocity. 



IV. EXCITON INTERCONVERSION PROCESSES 

Orthoexcitons in CU2O are observed to convert rapidly to paraexcitons [|l2j through an acoustic phonon- assisted 
process, depicted in Figs. 5(a,b), dependent on temperature but independent of the gas density at low densities [pT[. 
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At the zone center the orthoexcitons have T^ 5 symmetry and the paraexcitons Ig . Therefore, conversion requires a 
^25 = Tj 5 phonon. Since the acoustic phonons at the center of the zone in CU2O have rj~ 5 symmetry, the matrix 
element for the conversion process is strongly suppressed. Away from the zone center the matrix element, calculated 
in the k ■ p approximation, has the form |24j] , 

\M C \ 2 ~ k 2 , (17) 

where k a is the orthoexciton wavevector. In this approximation the symmetry of the orthoexciton wavefunction is 
approximately that at the zone center, plus a small component with rj~ 5 symmetry of amplitude cx k a . We begin by 
calculating the down-conversion rate (Figs. 5(a,b)) when the system is not in the condensed region. We extend the 
calculation to condensed excitons below. 
The down-conversion rate, D, is given by 



- tn^ c ? i M ci 2 K(i + + + A - H-t - 

k,q 

+n£/_ ? (l + nl_ ? )5{e % + A - e^+Huf)}, (18) 

where N 0texc is the number of non-condensed orthoexcitons. We extract the matrix element here from experimental 
data for down-conversion in the non-degenerate regime [Q, for which the classical limit of Eq. ( |lS| ) applies. We find 
for low lattice temperatures that in this limit, D ss 0.014T 3 / 2 ns -1 with T measured in Kelvin. To determine the rate 
of energy change rate due to down-conversion we multiply the first term in the summand of Eq. (Es|) by (A — hu^) 
and the second by (A + Jiu>$) . 

The up-conversion rate, U , for the processes shown in Figs. 5(c,d), is similarly, 



U 



1 / dN, 



+n l+^-^( 1 + + A " e s- f - ^?)>. ( 19 ) 

where N Ptexc is the number of non-condensed paraexcitons. Since the ortho-para energy splitting AE — 12 meV is 
^> fcgT for temperatures of interest, U is much smaller than D. To determine the rate of energy change due to this 
process we multiply the first term of the summand of Eq. (|l9|) by — (A + Tiuj^) and the second by (—A + Tioj^). 



We turn now to the condensed excitons. Equation (17) implies that down-conversion of condensed orthoexcitons 



(k Q = 0) is suppressed by the symmetry of the excitons and the phonons at the center of the zone of CU2O. The 
up-conversion of condensed paraexcitons is not possible either, because energy and momentum cannot be conserved 
for the parameters of CU2O. Thus only excitons in excited states participate in this mechanism. The conversion 
process primarily influences the orthoexcitons, since the up-conversion rate of paraexcitons is small. Condensation of 
orthoexcitons would in fact help them to maintain their high density and lower the heating by down-conversion. 

Figure 6 shows the calculated rates of the down and up-conversion processes as a function of temperature at fixed 
a. As expected both the down-conversion and the up-conversion rates increase with exciton temperature. On the 
other hand, both rates are suppressed as the degeneracy increases. 



V. AUGER PROCESS 

The non-radiative direct and phonon-assisted Auger decay processes, illustrated in Figs. 7,8, are the most important 
mechanisms for loss of excitons at high densities. Studies of this mechanism have been carried out on excitons in 
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CU2O P]l3|], and in electron-hole plasmas [[25j]26j where an electron recombines and excites either another electron 
high in the conduction band, or a hole deep in the valence band. Detailed calculations of the Auger rates are given 
in Ref. ||; we summarize the results here. 

The band structure of the material plays an essential role in the Auger annihilation of excitons. Because the 
conduction and the valence bands have the same (even) parity |2^J, the rate of the direct Auger process (Fig. 7) 
without Umklapp is negligible, and the Umklapp terms are small. Although the phonon-assisted Auger process 
(Fig. 8) is suppressed by the additional matrix element for phonon emission, the process is dominated by a large 
non-Umklapp term when the phonon has negative parity. The process is in fact faster than the direct. The rate 
of phonon-mediated para-para Auger annihilations is negligible compared with the rates for ortho-ortho and ortho- 
para; the reason is that the intermediate negative-parity band which enters the process is very close in energy to the 
conduction band for the orthoexciton recombination vertex entering ortho-ortho or ortho-para collisions, but it is a 
deep valence band for the paraexciton vertex entering para-para collisions. 
The Auger decay times appearing in Eqs. and (||) are given by || 

(20) 
(21) 

for the phonon-assisted process, where the constant C is of order 0.4 ns _1 , when the exciton densities, n^, are measured 
in units of 10 18 cm -3 . The Auger decay rate depends only on the density of excitons for low lattice temperatures. 
Since the Auger matrix element is constant when the thermal energy of excitons (a few meV) is much less than the 
gap energy E g (=2.17 eV in CU2O), the decay rate depends on neither the exciton temperature nor the statistics. 
Auger collisions are not influenced by the thermal motion of the excitons, and hence, the decay rate is the same for 
degenerate and classical gases of the same density. 

The Auger process produces hot carriers within the exciton gas - ionized electrons and holes - excited by an energy 
of order the energy gap. The hot carriers primarily lose this energy by emission of a cascade of phonons (Fig. 9(a)) - 
and to a lesser extent by scattering with excitons (Fig. 9(b)) - and eventually reform excitons p8fl . To calculate the 
fraction of Auger energy which heats the exciton gas we split the problem into two parts, before and after the exciton 
is formed. We shall see that only a small fraction of the hot-carrier energy goes into heating the gas, because the 
carriers are quickly cooled by rapid Frohlich emission of optical phonons. In contrast the hot excitons have a slower 
phonon emission rate, and their kinetic energy is more readily shared with the thermalized excitons. 

The cooling of the carriers by optical phonon emission (Fig. 9(a)) is given by 



/ dE \ 2 

opt 



V dt 



'i 



At high energies the Frohlich polar interaction |lq| is the dominant process for the emission and absorption of free 
carriers by optical phonons. If there is only one optical phonon branch, the matrix element for the Frohlich interaction 
is, 

\M,?^{^-^. (23) 

Here £oo = 6.46 is the high frequency and eo(= 7.11) is the static dielectric constant; these two dielectric constants 
differ only slightly as a consequence of CU2O being almost non-polar plfl . Since CU2O has a large number of optical 
phonon modes, Eq. (^3|) can only be applied approximately; for the lowest optical phonon of energy 11.4 meV, it 
leads to a rate ~ 50 eV/ns for the highest energy of the carriers (« 1 eV); the net rate should be even higher. 

The rate of energy change of a hot carrier of momentum k and energy E? due to collisions with thermalized excitons 
of momentum p and energy ep (Fig. 9(b)) is given by 
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dE \ 2 

-Qf ) =-yE \W I V(l + W)(^ S - £jf_ f )«(e^+ - (24) 
where M| e is the carrier-exciton matrix element. If we assume a momentum transfer q ^> p, then 



(25) 



9* 

As we derive in the Appendix, the matrix element M£ e is 

27 rey (A e -A fe ) 
9 eo^ (1 + (ga/4) 2 ) 

where^ A e = m e /m = 0.58 and = m/j/m = 0.42. This process gives a cooling rate ~ 1 eV/ns, some 2% of the 
lowest optical phonon cooling rate, Eq. (p2|). Furthermore, the ratio of the rates of carrier-exciton to phonon cooling 
decreases with decreasing carrier energy. 

When the hot carriers lose sufficient energy they reform hot excitons which heat the exciton gas by scattering with 
thermal excitons. It is this heating of the exciton gas which stands in the way of the excitons crossing the phase 
boundary into the condensed region. Cooling of hot excitons by phonon emission is competitive, but does not heat 
the exciton gas. Heating of the excitons per Auger pair recombination is at most of order the exciton binding energy. 
To calculate the kinetic energy transfered to the exciton gas by this Auger annihilation process, we note that the 
newly formed excitons scatter against the thermalized excitons (Fig. 9(c)) and cool at a rate 

de h \ e h e h 

—^=-- n v rd o[v rd ) = - — , (27) 

where et is the kinetic energy of a hot exciton, a(v re i) is the exciton-exciton scattering cross section, v re i the relative 
velocity of the excitons, and n the thermalized exciton density. At low energy the cross-section is cr(v re i) — 4ir(2a sc ) 2 , 
with a single-exciton radius for exciton-exciton collisions, a sc . We use a sc — 3.0A which approximates the numerically 
calculated interaction potential p9| quite well. At high energy the cross section calculated in Born approximation 
taking into account Coulomb interactions between the electrons and holes in the different excitons is 

33 / e 2 \ 2 

a(v rel ) = —na 2 — . (28) 

35 \riVreieoJ 

A simple interpolation between these limits is 

4^(2a sc ) 2 
1 + (140/33) {k re ia sc y 

where k re i = (mv re i/h). In the limit that the hot exciton energy is much larger than that of the thermalized excitons, 
Eqs. d^||) give 



dtj exc l \aJ {Tilo,) 2 \ma 2 J {1 + {280/33)ma 2 sc e h /h 2 ) 



The prefactor, huf is « 198 eV/ns for the lowest energy optical phonon of 11.4 meV. 

Hot excitons also cool by phonon emission (Fig. 9(d)) without heating the exciton gas. The rate of cooling of a hot 
exciton of momentum k and energy by optical phonon emission is given by 



2 As mentioned earlier the exciton mass is not equal to the sum of the effective electron and hole masses; we assume, however, 
that the ratios m e /m and rrih/rn are the same as those following from the assumption m = m e + nih- 
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= - ¥E W(e S - *s_ ? - (31) 

at / opt " f 

with squared matrix element ( |l3|) . Emission of optical phonons via deformation potential coupling cools a hot exciton 
of kinetic energy eh > fiioi at a rate 



V2 m 2 m „-^. 1/2 / fi 2 X 1/2 

(32) 



^ /opt i 1 J m ion \ aJ 2 1 /2 7r (7j a;i )3/2 y ma 

For hot excitons, all 15 optical phonon modes need to be taken into account in the sum. The measurements of Snoke 
et al. probe only the three lowest optical phonons T 25 , T^ 2 and rj~ 5 , of energies 11.4, 13.8 and 18.7 meV, respectively. 
For the rj~ 2 optical phonon of 13.8 meV, for example, 

de h 
dt 




-0.16 — 1 eV/ns. (33) 



In addition, excitons of energy higher than Ihioi can cool by a fast parity-conserving two optical-phonon emission 
process [fl8f . For the simultaneous emission of two T^ 2 phonons of 27.6 meV, Snoke et al. have placed experimental 
limits on the cooling rate which indicate that it is at least ten times larger than the single phonon rate. They have 
also found that the two-r^ phonon cooling mechanism is too small to be observed. The complete treatment of optical 
phonon cooling including all phonon processes is beyond the scope of our present study and indeed requires much 
more experimental input. Nevertheless, we can use the experimentally-derived exciton cooling rates of Snoke et al. 
involving the three lowest optical phonon branches as a gauge of the total cooling rate for all optical phonon processes. 
We introduce the cooling parameter 7 (> 1) such that 



de h 
dt 



= 7 

opt 



de h \ ( den 

dt L- \ dt 



2rr 



(34) 



We shall see that a fit of the short-pulse experimental results yields 7 « 4.5. 

At energies lower than the lowest optical phonon energy only acoustic phonons can be emitted; the cooling rate for 
coupling to longitudinal acoustic phonons calculated with a deformation potential coupling is 

__^J!L(!£)'_2^fJ?,) ,/, S&M, (35) 

dt ) la m ion \aJ Sirihtutfe]/ 2 \ ma e J (1 + z ) 3 

where zo = 2 [(2me/ l a 2 ) 1 / 2 — (mi^a)] /%. In general the exciton-phonon interaction is not as effective at cooling 
excitons as it is as for free carriers. For the polar interaction this is due to the electrical neutrality of the exciton and 
for the deformation potential scattering the matrix element of the process involves the difference of the deformation 
potentials for the conduction and the valence bands. 

The total cooling rate of a hot exciton due to phonon emission and scattering with thermalized excitons is 



deh _ ( de h \ ( den 
dt \ dt) ph \ dt 



(36) 



where (deh/dt) p h is the total phonon cooling rate of hot excitons, the sum of Eqs. (Q) and (p5|). The heat transfered 
to the exciton gas, £U, per Auger annihilation process is then 

where to is the time it takes for the hot exciton to fall to the mean kinetic energy of a thermalized exciton. This time 
must be found numerically by solving for the time evolution of the kinetic energy of the hot exciton. 



10 



We note that Auger heating depends strongly on the exciton density; it depends only weakly on the exciton 
temperature because the higher the kinetic energy of the exciton gas, the less energy is transfered from hot to cold 
excitons. Auger heating does not depend on the degeneracy of the gas. In the limit that the density of the thermalized 
excitons goes to zero, the Auger heating goes to zero. In the limit that the exciton density is very large, the Auger 
heating is equal to the binding energy Ef,, since all the energy goes to the excitons. The ratio Ea/E/, as a function of 
the density of the thermalized excitons is shown in Fig. 10, assuming that the exciton gas temperature is kept fixed 
at T = 30K. 

The rate of energy change of the exciton gas due to Auger heating is 

1 fdE\ E A 



N\dt) A ~ T f (38) 

where {rfy 1 = -{l/N){dN/dt) A = {C/A)n is the total Auger scattering rate for the excitons, Eqs. ( pp| , ^l|) . Since 
the laser creates ionized electrons and holes with an excess kinetic energy of about 250 meV per pair, an energy larger 
than Eb, the heating of the exciton gas from hot carriers generated by the laser is the same as the Auger heating 
Ea- The rate at which the laser heats the excitons is effectively Ea times the laser production rate, and for a given 
density and temperature of the exciton gas the total heating rate due to Auger decay and the laser production is E A 
times the sum of the Auger rate and the laser generation rate. 



VI. RESULTS 

In the previous sections we have derived from experimental measurements and theoretical considerations numerical 
values of the parameters entering the rates of energy change of the exciton gas due to phonon cooling, the rates 
of ortho-para interconversion and energy exchange, and the Auger decay rate. The processes of Auger and laser 
heating involve the uncertain cooling parameter 7, Eq. j3^). We now solve the kinetic equations for N a (t), N p (t) 
and T(t) numerically, assuming reasonable initial conditions for these quantities at t = 0. We choose 7 to optimize 
agreement between our results and the experimental data, for the case of short-pulse excitation, (subsection A; Fig. 
11) and keep this parameter fixed in the rest of our calculation. We find that 7 w 4.5, a reasonable number as an 
average of the cooling rate of a hot exciton due to all the optical phonons of CU2O (including single and multi-phonon 
processes) in terms of the r^ 2 plus 2T^ 2 optical-phonon cooling rate. 

In the present calculations of the dynamics we assume that the orthoexcitons are triply degenerate and, hence, we 
compare our orthoexciton calculations with the zero-stress data. No experimental determination of the paraexciton 
density has been made in the zero-stress case. We find, however, that predictions for the paraexcitons at zero stress 
agree remarkably well with the observed paraexciton behavior under stress. This is not surprising as the paraexciton 
multiplicity does not change with applied stress. 

The numerical simulations together with the argument, given in the Introduction, that the orthoexcitons should 
approach an adiabat just above the condensation line, put us in a position to identify the crucial factors that lead 
to the quantum saturation of orthoexcitons, while allowing condensation of paraxcitons. At the outset, we note that 
no single process we consider poses an obviously insurmountable barrier to Bose-Einstein condensation. To show in 
detail the physical effects of the various processes we take the approach of first computing the behavior of the gas 
with all the known factors in the equations, and then dissecting the problem by removing each of the processes from 
the calculation in turn. 



A. "Short-pulse" excitation 

A mode-locked, cavity-dumped Argon-ion laser produces nanojoule pulses with about a 100 ps length. For the 
photoluminescence time resolution of about 100 to 300 ps, this excitation pulse is effectively a delta-function in time, 
and the evolution of the system is observed without creation by the laser of further particles. The excitons, whose 
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energy relaxation and decay we are dealing with, form on a timescale shorter than the detection limits; as we shall see, 
most of the action occurs within a few nanoseconds. We compute the behavior of both orthoexcitons and paraexcitons 
in this "short-pulse" case, even though the radiative efficiency of the paraexcitons is a factor 500 less than that of the 
(higher energy) orthoexcitons, making difficult to observe the paraexcitons during this short time period. 

Figure 11 shows the result of our simulation with initial conditions, n = 25 x 10 18 cm -3 and n p = 8.33 x 10 18 cm -3 
at a temperature of 60K, and Ti fixed at 8K. We assume in this and the long-pulse case that the laser produces three 
orthoexcitons for each paraexciton, as discussed earlier. For comparison to experiment in this figure and those that 
follow, we plot Snoke's orthoexciton data (diamonds) for the long-pulse case between n a = 10 18 cm -3 and 2 x 10 19 
cm~ 3 . The solid and dashed curves represent the theoretical orthoexciton and paraexciton trajectories, respectively, 
over the time t = to 2 ns. As time progresses the excitons cool. The dotted straight lines are the Bose-Einstein 
condensation lines, the upper for paraexcitons (g = 1) and the lower for orthoexcitons (g = 3). 

A principal result of our study is that, as in experiment, the orthoexcitons move closely parallel to the phase 
boundary without crossing it, following an adiabat, as a result of the balance between the phonon cooling and the 
Auger heating, Eq. ([!]). We present the details of this argument in subsection C. When the orthoexciton density falls 
sufficiently low (< 10 18 cm -3 ), and their temperature approaches that of the lattice, the orthoexcitons move away from 
the condensation line; their density continues to decrease due to ortho-to-para conversion and Auger recombination. 

The calculations further predict that the paraexcitons quickly cool below their phase boundary and condense. There 
is little decay in the density of the paraexcitons over the 2 ns time interval; in fact, their density initially increases 
due to the down-conversion of orthoexcitons. Once the orthoexcitons are gone, the paraexcitons decay on a scale of 
their recombination time, which is longer than the chosen time interval. Note that the two components cover the 
same temperature range because we have assumed, as throughout, that they are in good thermal contact with each 
other and thus have the same temperature. 

Now let us examine how changing specific quantities in our calculation changes the results, starting with the lattice 
temperature. We expect that at lower lattice temperatures, the orthoexcitons cool more efficiently. We see from 
Fig. 12(a) that this is true to a small extent. With a lattice temperature of 4K, the orthoexcitons still follow the 
same trajectory at high density - implying that the heating is associated with other processes - but as the gas 
density reduces below about 3 x 10 18 cm~ 3 a brief crossing of the phase boundary occurs. Unfortunately, there is 
no direct experimental control over the local lattice temperature: the bath temperature is already 2K and the rise 
in lattice temperature is indirectly caused by the phonon emission of photoexcited carriers. Matters could be worse: 
we see in the figure that a lattice temperature of 18K makes Bose-Einstein condensation highly improbable for the 
orthoexcitons. 

The assumption of quantum statistics in the exciton-phonon cooling processes has a marked effect on the quantum 
saturation, as shown in Fig. 12(b). Taking classical statistics (hypothetically) in the cooling process, we see that 
the orthoexciton gas does not follow a constant-a trajectory at any stage. The quantum effects shown in Fig. 4 - 
reductions in the cooling rates at high quantum degeneracy - apparently have a deleterious effect on the condensation. 

Now let us restore the quantum statistics and observe the effect turning off ortho-to-para conversion. As we see 
in Fig. 13 the paraexcitons, while still predicted to condense, do not show an increase in density with time, because 
one source of paraexcitons has been removed. Still, Auger recombination diminishes the density of both types of 
particles. The orthoexcitons, however, come considerably closer to their phase boundary because one of their decay 
mechanisms has been deleted. Yet, no condensation is predicted for the orthoexcitons. Ortho-to-para conversion is 
not the dominant culprit in preventing their condensation. 

Next let us see how the Auger recombination process affects the system. In the calculation shown in Fig. 14(a) 
we assume that the Auger rate (and, therefore, the Auger heating) is reduced by 50%. This decrease enhances the 
density of both species at a given time and lowers the temperature of the exciton gas. We see that the orthoexcitons 
are barely able to condense in this case. To carry this effect to the extreme we remove the Auger recombination 
completely, with the result shown in Fig. 14(b). In this case the total number of excitons is nearly conserved on the 
2 ns timescale because the radiative lifetimes of the excitons are much longer. Ortho-to-para conversion causes the 



12 



density changes seen in the figure. The condensation of both the orthoexcitons and paraexcitons is now unavoidable. 
Actually there are two reasons for this: first, the decay of both species is reduced, and, second, there is no Auger 
heating, which depends strongly on the exciton density and pushes the excitons away from the phase boundary. 

Optimistically, we can say that Nature has been kind. Were the Auger rate 50% larger, it is unlikely that even the 
quantum saturation effect - which, as we have seen here, indeed senses the quantum statistics - would be observable. 

B. "Long-pulse" excitation 

If the mode-locker is removed from the laser, the cavity-dumped mode provides 10 ns-long pulses with about an 
order of magnitude more energy. Now we must incorporate the generation rate of the laser in the numerical calculation. 
We assume that the initial temperature of the exciton gas is equal to the lattice temperature, which we take to be 
6K. To start the gas off in the classical regime, we choose initial densities of 0.1 x 10 18 cm -3 for the orthoexcitons 
and 0.033 x 10 18 cm~ 3 for the paraexcitons. We further assume that the temporal laser profile is Gaussian and as 
before that three orthoexcitons are produced for each paraexciton. 

The principal result of our calculation is shown in Fig. 15. Now the orthoexcitons initially move up the phase 
boundary, towards higher density and temperature, and decay closely along the same path. This is precisely the 
remarkable quantum-saturation behavior which we hoped to reproduce theoretically and understand. In addition, 
we see that the paraexcitons heat up with the orthoexcitons, but they are able to cross their phase boundary at an 
early stage, implying that the gas has a condensed fraction. At the peak of the laser pulse, both orthoexciton and 
paraexciton densities begin to decrease. The paraexcitons retain a high density as they cool and eventually reach a 
condensed fraction of over 90%! The numerical simulation is remarkably similar to the data of Lin and Wolfe in 
their uniaxially stressed sample. As previously stated, we do not expect major changes in the calculated behavior for 
the stressed and unstressed cases. 

For clarity, we reproduce in Figs. 16 the results of Fig. 15, showing separately the build-up and decay processes. 
In our calculation, a pulse with a full-width at half-maximum of 3 ns is used, as shown in the insets of Figs. 16. 
Figure 16(a) plots the (upwardly rising) gas trajectories for the first half of the laser pulse, and Fig. 16(b) shows the 
(downward falling) results for the second half of the laser pulse. We note that the initial conditions of the "decay 
phase" differ from those of the short-pulse experiments: here, the "initial" paraexciton density is considerably higher 
than the orthoexciton density, because the ortho-to-para conversion process has had the entire first half of the laser 
pulse to pump up the paraexciton density. 

Again, we test the effect of quantum statistics in the phonon cooling processes. Figure 17 shows that in the case 
of "classical" phonon-cooling, both orthoexcitons and paraexcitons are predicted to condense. Thus, the quantum 
statistics play a crucial role in keeping the orthoexcitons away from the condensed region. 

In Fig. 18 we show the result of our calculation if we neglect the ortho-para interconversion mechanism. As seen 
from this figure, the orthoexcitons condense, but still tend to move along the condensation line. A major difference 
between this graph and the data shown in Fig. 2(b) is the faster decay rate here of the paraexcitons at late times. 
The reason for this is the lack of the ortho-to-para conversion mechanism, which balances the loss of paraexcitons due 
to the Auger process. 

In Fig. 19 we show the result of our calculation if the Auger process (annihilation and heating) is neglected. Both 
species condense for the same reasons as those given for the short-pulse excitation, in Fig. 14(b). The effect of 
the Auger heating itself (with unmodified Auger annihilation) is shown in Figs. 20 (a, b), where we assume that the 
scattering radius a sc for carrier-carrier scattering is ten times smaller or larger, respectively, than assumed in Fig. 15. 
This scattering length influences the kinetic energy transfcrcd from the Auger-ionized carriers to the thermalized 
excitons. In the first case both species condense because of the low Auger heating. In the second case, however, the 
behavior of excitons hardly changes from that of the initial calculation, Fig. 15. The reason for this is that, according 
to Eq. (^9|) the cross section for the scattering is independent of the exciton radius at large radius, and in this limit 
the rate of energy transfer is dominated by the scattering process. 
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Finally in Fig. 21 we investigate the result of changing the Auger decay constant C and the acoustic phonon 
deformation potential D a . The constant C determines the Auger decay rate and consequently the Auger heating rate. 
Also D a determines the dominant phonon cooling rate, which is the acoustic phonon mechanism. For values of C and 
D a inside the shaded region, the orthoexcitons move along and closely to the condensation line, i.e., between the two 
critical lines for condensation of the ortho and the paraexcitons. The region corresponds to < a < 0.45. For values 
of C and D a above the shaded region, the orthoexcitons condense because of the more effective phonon cooling, the 
low Auger decay rate and the low Auger heating. Below the shaded region the orthoexcitons move away from the 
degenerate regime, towards the classical regime, because of the low phonon-cooling rate, the high Auger decay rate 
and the high Auger heating. 



C. The entropy of orthoexcitons 

In the Introduction we briefly outlined the argument that the balance between the rates of phonon cooling and 
Auger heating causes the orthoexcitons to approach a critical adiabat, a line of constant entropy per particle of the 
orthoexcitons, s 0l lying slightly above the condensation phase boundary. We now spell out the calculation in detail. 
The rate of change of s a is given by 

dsg = 1 dE Ts + Ho dN P dV 

dt N dt N dt N dt ' 1 ' 

where E Q is the total energy and P the pressure of the orthoexciton gas. 

The dominant contributions to the right side of Eq. (|39| ) are the phonon cooling and Auger heating, which enter 
the first term. The second term is negligible, since particle numbers change relatively slowly. For short pulses both 
the ortho-to-para conversion and the Auger decay rates are proportional to n ~ y 3 / 2 along the phase boundary; the 



total decay rate is w 1 ns" 1 for (orthoexciton) density 10 18 cm -3 |13|. For the case of long pulses the net production 
rate (laser production minus decay processes) is even smaller. For small values of the chemical potential, therefore, 
the second term will be on the order of the kinetic energy times the decay rate, which is more than an order of 
magnitude smaller than the phonon cooling and the Auger heating rates. In addition, the final term is negligible for 
small expansion of the gas, which we have assumed. 

The slow rate of change of particle numbers also allows us to relate the first term on the right side of Eq. ( |39| ) to 
the total change of the energy of the excitons (ortho plus para) by 

1 dE N C (T) 1 dE 



N dt N C (T) + C P {T) N dt ' 

where Cj = (dEi/dT)isr it v is the heat capacity of the i th component at constant volume and number of particles. For 
the excitons close to the Bose-Einstein condensation phase boundary, C p /C ~ N p /N 0l so that (l/N )(dE /dt) « 
(l/N)(dE/dt). 

The acoustic-phonon cooling rate per exciton, for excitons with a = 0.15, is w — 0.56T 3 / 2 (1 — T(/T) meV/ns [cf. 
Fig. 4a]. Since we are interested in T 3> Te, we neglect the term T&/T, 

In the case of short-pulse excitation, the paraexciton density is approximately constant, « 10 19 cm -3 , due to a 
balance between the ortho-to-para conversion process and the Auger loss of paraexcitons. Since the Auger heating of 
orthoexcitons per particle is a function of the total exciton density (xL 10 19 cm" 3 ), we assume that Ea ~ 0.77E o [cf. 
Fig. (10)] and neglect its weak density dependence. Then the Auger heating rate per particle varies approximately 
as 11.8n c meV/ns with n a measured in units of 10 18 cm -3 . For the case of long-pulse excitation, the Auger heating 
is a bit smaller, closer to Ea ~ 0.64-E&, since the total density is smaller in this case, but still the density dependence 
of Ea is weak for n > 5 x 10 18 cm" 3 . 

With these approximations Eq. (^9|) takes the form 

T^«-aT 3 / 2 + 6n , (41) 
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where, for T in Kelvin and n Q in units of 10 18 cm -3 , a w 0.56 meV/ns, and b sa 11.8 meV/ns for short pulses and 
w 9.8 meV/ns for long pulses. Note that the right side of Eq. ( fi"l| ) vanishes on the adiabat for which, in the same 
units, n ss 0.048T 3 / 2 for short pulses, and n D w 0.058T 3 / 2 for long pulses, and that this equation predicts that the 
orthoexcitons approach this adiabat as a stable fixed point. The adiabat of stability, as calculated approximately 
here, is shown as a dash-dot line in Fig. 11 for short pulse-excitation and Fig. 15 for long-pulse excitation. 

To calculate the characteristic time, r», for the orthoexcitons to approach the adiabat of stability, we regard Eq. 
( [i"l| ) as an equation for a. Near this adiabat, on which a = a*, 

da _ a OLjf 

dt~ n ' { ] 

where 

__T(dsA ^ 5C(5/2) k B 

b \dn ) T ~ 2C(3/2)aTV2' ^ 

the latter form holds near the phase boundary. Numerically t* ~ 0.20T -1 / 2 ns, with T measured in Kelvin. In the 
case of short-pulse excitation, for T = 60K, t* w 25 ps. The present data is consistent with the predicted approach 
to the adiabat of stability, but since the photolumincscence time resolution in the experiments, from 100 to 300 ps, is 
much larger than r*, the approach cannot be verified in detail. 



VII. SUMMARY 

We have shown that the observed quantum saturation of orthoexcitons is caused by a competition between the 
Auger process and the cooling of the excitons by the phonons. When the orthoexcitons approach the phase boundary, 
these two mechanisms act together to keep the orthoexcitons on a particular adiabat, parallel, but at slightly higher 
temperature than the phase boundary. The cooling is affected by the proximity of the excitons to the phase boundary, 
but the Auger process is not significantly affected by the degeneracy of excitons. We find that there is a fairly wide 
range of numbers that can be inserted for the various parameters and still account for the observed quantum saturation 
effect. Reasonable changes in the parameter values cause the trajectories to shift relative to the phase boundary but 
do not qualitatively change the results. Also, the assumed laser profile and width do not change the results of our 
calculation appreciably. We have, however, made the assumption of a constant gas volume over the short time intervals 
considered (1-20 ns). 

Our key finding is that, while the orthoexcitons encounter a formidable barrier to condensation, the paraexcitons 
do indeed cross the phase boundary and condense. Specific features of the band structure of CU2O underlie this 
condensation. Earlier work has noted the simple parabolic bands and the forbidden direct gap. Our present study 
points out that paraexciton-paraexciton Auger collisions are very slow for both the direct and the phonon-assisted 
Auger processes. Details will be reported in Ref. Finally we remark that this analysis presents no fundamental 
reason why the orthoexcitons cannot condense; however their multiplicity, their faster Auger decay and their conversion 
to paraexcitons make their condensation much more difficult than for paraexcitons. 
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APPENDIX A: CALCULATION OF ELECTRON (OR HOLE)-EXCITON SCATTERING MATRIX 

ELEMENT 



We calculate here the matrix element M£ e for scattering of an electron or hole by an exciton (Fig. 9(b)). Consider 
an electron of momentum k and an exciton of momentum p scattering into k — q and p + q, respectively. The electron 
(or hole)-exciton interaction is: 

'!■ (Al) 



eo 



where f is the position of the electron (or hole) and f e and fh are the positions of the bound electron and hole, 
respectively. The wavefunctions of the initial and final states are 



rh), 



where 



^ = ^e^-% + ^f e -r ft ), 



%(r e - fh) = -l=e^ XA + XA ^ rel (re - fh) 



(A2) 



(A3) 



is the exciton wavefunction with center-of-mass momentum p, and <fr re i(f e — fh) is the wavefunction of the relative 
motion. The matrix clement of the interaction between the initial and final states is thus 
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[1 + (qa\ h /2)Y [1 + (gaA e /2) 
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(A4) 



The factor inside the parentheses comes again from the Fourier transform of the relative electron-hole wavefunction, 
which we assume to be the ls-state. In the limit A e w A/,, 



Mf 



27re 2 a 2 (A e - X h ) 
eoV (1 + (ga/4) 2 ) 3 



(A5) 



[1] T.J. Greytak and D. Kleppner, in New Trends in Atomic Physics, G. Grynberg and R. Stora eds. (North Holland Publ., 
Amsterdam, 1984). 

[2] Bose Einstein Condensation, A. Griffin, D.W. Snoke, and S. Stringari, eds. (Cambridge Univ. Press, Cambridge, 1995). 
[3] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, and E.A. Cornell, Science 269, 198 (1995); C.C. Bradley, 

C.A. Sackett, J.J. Tollett, and R.G. Hulet, Phys. Rev. Lett. 75, 1687 (1995); K.B. Davis, M.-O. Mewes, M.R. Andrews, 

N.J.van Druten, D.S. Durfee, D.M. Kurn, and W. Ketterle, Phys. Rev. Lett, (in press). 
[4] D. Hulin, A. Mysyrowicz, and C. Benoit a la Guillaume, Phys. Rev. Lett. 45, 1970 (1980); A. Mysyrowicz, D. Hulin, and 

C. Benoit a la Guillaume, J. Luminescence 24/25, 629 (1981). 
[5] D.W. Snoke, J.P. Wolfe, and A. Mysyrowicz, Phys. Rev. Lett. 59, 827 (1987). 

[6] D.W. Snoke, J.P. Wolfe, and A. Mysyrowicz, Phys. Rev. Lett. 64, 2543 (1990); Phys. Rev. B 41, 11171 (1990). 
[7] J.-L. Lin and J.P. Wolfe, Phys. Rev. Lett. 71, 1223 (1993). 

[8] V.A. Kiselev and A.G. Zhilich, Fiz. Tverd. Tela 13, 2398 (1971) [Sov. Phys.-Solid State 13, 2008 (1972)]; M.M. Denisov 
and V.P. Makarov, Phys. Stat. Sol. B 56, 9 (1973); V.P. Makarov, Zh. Eksper. Teor. Fiz. 54, 324 (1968) [Sov. Phys. JETP 
27, 173 (1968)]; G.E. Pikus and G.L. Bir, Zh. Eksper. Teor. Fiz. 60, 195 (195) [Sov. Phys. JETP 33, 108 (1971)]. 
[9] G.M. Kavoulakis and G. Baym, to be published. 
[10] We note, however, that even in the stressed case, the paraexciton spectrum is not completely explained. Its non-zero width 
and slight blue shift have been modeled by H. Shi, G. Verchaka and A. Griffin [Phys. Rev. B 50, 1119 (1994)] in terms 
of many-body interactions within the gas. Understanding this spectrum remains one of the outstanding problems of the 
quantum-degenerate exciton system. 



16 



[11] D.W. Snoke, D.P. Trauernicht, and J. P. Wolfe, Phys. Rev. B 41, 5266 (1990). 

[12] F.I. Kreingol'd and V.L. Makarov, Fiz. Tverd. Tela 15, 1307 (1973) [Sov. Phys.-Solid State 15, 890 (1973)]; J.S. Weiner, 

N. Caswell, P.Y. Yu, and A. Mysyrowicz, Sol. State Commun. 46, 105 (1989). 
[13] D.W. Snoke and J.P Wolfe, Phys. Rev. B 42, 7876 (1990); D.P. Trauernicht, J. P. Wolfe, and A. Mysyrowicz, Phys. Rev. 

B 34, 2561 (1986). 
[14] B. Link and G. Baym, Phys. Rev. Lett. 69, 2959 (1992). 
[15] J. Bardeen and W. Shockley, Phys. Rev. 80, 72 (1950). 

[16] K. Seeger, in Semiconductor Physics: an introduction, Springer series in solid-state sciences, v. 40 (Springer- Verlag, New 
York, 1991). 

[17] P.Y. Yu and Y.R. Shen, Phys. Rev. B 12, 1377 (1975). 

[18] D.W. Snoke, D. Braun, and M. Cardona, Phys. Rev. B 44, 2991 (1991). 

[19] M. O'Keefe, J. Chem. Phys. 39, 1789 (1963). 

[20] A.G. Zhilich, J. Halpern, and B.P. Zakharchenya, Phys. Rev. 188, 1294 (1969). 

[21] JW. Hodby, T.E. Jenkins, C. Schwab, H. Tamura, and D. Trivich, J. Phys. C 9, 1429 (1976), T. Itoh and S. Narita, J. 

Phys. Soc. Japan 39, 132 (1975). 
[22] C. Carabatos and B. Prevot, Phys. Stat. Sol. B 44, 701 (1971), K. Reimann and K. Syassen, Phys. Rev. B 39, 11113 

(1989). 

[23] D.P. Trauernicht and J.P. Wolfe, Phys. Rev. B 33, 8506 (1986). 

[24] E.O. Kane, in Semiconductors and Semimetals, v. 1, Physics of III-IV Compounds, R.K. Willardson and A.C. Beer, eds. 

(Academic Press, New York, 1966), p. 75. 
[25] P.L. Gourley and J.P. Wolfe, Phys. Rev. B 24, 5970 (1981); H. Yoshida, H. Saito, and S. Shionoya, Phys. Stat. Sol. B 104, 

331 (1981). 

[26] A.R. Beattie and P.T. Landsberg, Proc. Roy. Soc. A 249, 16 (1958); B.K. Ridley, in Quantum Processes in Semiconductors 
(Clarendon Press, Oxford, 1982); A. Haug, Sol. State Commun. 22, 537 (1977); A. Das and R. Al-Jishi, Phys. Lett. A 
141, 186 (1989); A. Haug, Sol. State Electr. 21, 128 (1978). 

[27] J.P. Dahl and A.C. Switendick, J. Phys. Chem. Solids 27, 931 (1965); R. J. Elliott, Phys. Rev. 124, 340 (1961); L. 
Kleinman and K. Mednick, Phys. Rev. B 21, 1549 (1980). 

[28] S.E. Esipov, M. Msall, and J.P. Wolfe, Phys. Rev. B 47, 13330 (1993). 

[29] E. Hanamura and H. Haug, Phys. Rep. C 33, 209 (1977). 



FIG. 1. (a)-(c) The LO phonon- assisted recombination spectrum of orthoexcitons for three different times following 
a 100 ps laser pulse. The data are fitted by Bose-Einstein distributions (dashed lines). Temperatures and densities 
extracted from the fits to the photoluminescence spectra, (d) Open circles show the density of orthoexcitons at various 
times calculated from experimentally determined \i and T, as a function of temperature, for short laser pulses. Black 
dots are the results for orthoexcitons in the long-pulse (10 ns) case. In both experiments the orthoexcitons move 
along the critical line for condensation, without condensing. All data is from Snoke et al. Refs. |5||| and ]l3| , 

FIG. 2. (a) Data from lightly stressed crystals with long-pulse (10 ns) excitation Q. The laser profile (triangles), 
the number of orthoexcitons in the lowest orthoexciton level (open circles), and the number of paraexcitons (black 
dots) as function of time. The paraexcitons show a significantly smaller decay rate, (b) Corresponding trajectories 
for orthoexcitons (open circles) and paraexcitons (black dots) in the density-temperature plane. The straight line is 
the condensation phase boundary, which is identical for paraexcitons and orthoexcitons in the stressed case. Note 
that the paraexcitons are in the condensed region at times later than 8 ns. 

FIG. 3. Exciton-phonon interaction processes for non-condensed (a,b) and condensed (c,d) excitons. 

FIG. 4. Acoustic and optical phonon cooling rates per exciton as a function of a for fixed exciton gas temperature 
30K and lattice temperature 6K. Increased degeneracy (smaller a) lowers the cooling rate. 

FIG. 5. Orthoexciton-paraexciton (a,b) down-conversion and paraexciton-orthoexciton up-conversion (c,d) for 
acoustic phonon-mediated mechanisms. 

FIG. 6. The down-conversion and up-conversion rates are shown for a = a p = and a = a p = 0.1 as a function 
of the temperature of the exciton gas. 
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FIG. 7. Direct Auger non-radiative annihilation processes. Time progresses from left to right. The initial state 
contains two excitons of momenta K and P, and the final state contains an ionized electron and hole with momenta 
k e and fc/j, respectively. The dashed line denotes the Coulomb interaction. 

FIG. 8. Phonon-assisted Auger non-radiative annihilation processes. The wiggly line denotes a phonon of momen- 
tum Q. 

FIG. 9. Energy transfer processes among hot electrons and holes (single lines), excitons (double lines), and phonons 
(wiggly lines), (a) and (b) describe (Auger or laser generated) electron (hole)-phonon and electron (hole)-exciton 
scattering processes, respectively; (c) hot exciton-thermalized exciton scattering; (d) hot exciton-phonon scattering. 

FIG. 10. The ratio of the Auger heating over the binding energy Ef, versus the density of the thermalized 
excitons, at fixed exciton temperature, 30K. 

FIG. 11. Numerical solution of Eqs. (§^) in the density-temperature plane for the orthoexcitons (solid line) and 
paraexcitons (dashed line), for a short laser pulse. The diamonds are (long-pulse) experimental data for orthoexcitons. 
The two straight lines are Bose-Einstein condensation boundaries for orthoexcitons (lower) and paraexcitons (higher) 
The dash-dot line shows the adabiat of stability (see Eq. (|4l])) with only acoustic-phonon cooling and Auger heating 
taken into account. 

FIG. 12. (a) Calculated behavior of orthoexcitons for three different lattice temperatures, 18K (higher), 8K (middle) 
and 4K (lower), assuming a short laser pulse, (b) Behavior of orthoexcitons assuming classical and degenerate statistics 
for the exciton-phonon interaction. 

FIG. 13. Short laser-pulse excitation. No ortho-para interconversion process. 

FIG. 14. Short laser-pulse excitation, (a) Auger rate reduced by 50%; (b) Auger process turned off. 

FIG. 15. The behavior of ortho and paraexcitons in the density-temperature plane, for a long (10 ns) laser pulse. 
The dash-dot line shows the adabiat of stability (see Eq. (|4l])) with only acoustic-phonon cooling and Auger heating 
taken into account. 

FIG. 16. The build-up (a) and decay (b) phase of Fig. 15. 

FIG. 17. Behavior of ortho and paraexcitons assuming classical statistics for the exciton-phonon interaction, for 
long laser-pulse excitation. 

FIG. 18. Behavior of ortho and paraexcitons in the density-temperature plane, ignoring the interconversion process, 
for long laser-pulse excitation. 

FIG. 19. Behavior of ortho and paraexcitons in the density-temperature plane, ignoring the Auger (heating and 
decay) process, for long laser-pulse excitation. 

FIG. 20. Behavior of ortho and paraexcitons in the density-temperature plane, for long laser pulse excitation, 
assuming that the scattering radius for exciton-exciton collisions (a) reduced by a factor of ten, and (b) ten times 
larger than in Fig. 15. 

FIG. 21. The shaded region shows the range of the values of the Auger heating parameter C and the magnitude of 
the acoustic-phonon deformation potential D a which gives the observed quantum saturation of the orthoexcitons. 
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